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Abstract 

The tail problem for the propagation of a scalar field is considered in a cosmo- 
logical background, taking a Robertson- Walker spacetime as a specific example. 
The explicit radial dependence of the general solution of the Klein-Gordon equa- 
tion with nonminimal coupling is derived, and the inapplicability of the standard 
calculation of the reflection and transmission coefficients to the study of scattering 
of waves by the cosmological curvature is discussed. 
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1 Introduction 



The Klein-Gordon equation^ 

g ab V a V b ^ -£R$ -m 2 $ = , (1.1) 

where £ is a numerical constant, R is the curvature scalar, and m is the mass of the 
scalar field describes the simplest examples of relativistic wave propagation in curved 
spacetime and has consequently received much attention, especially in relation to the 
study of quantum processes 0. Nevertheless, some interesting features have not yet 
been fully investigated: In particular, the physical consequences of the fact that the 
solutions of Eq. (|1.1|) , even in the case in which £ = m = 0, generally present tails 
p|- ||10|| (i.e., propagation in the interior of the light cone), has been considered only in 
the context of radiation emitted by compact objects fl^-[|i~4|j. 

The formation of tails when £ = m = can be understood as an effect of the 
backscattering of curvature on the waves; this is schematically represented in Fig. 1, 
which represents a portion of spacetime in null coordinates u, v. The curvature is 
nonzero only in the shaded region, and scatters waves analogously to what happens 
under the action of a potential. The ray 1, which does not enter the curved region, 
remains unperturbed, whereas the ray 2 can be decomposed into a transmitted (3) and 
a reflected (4) part. Nevertheless we emphasize that such a representation is purely 
pictorial and does not properly describe many features of the real phenomenon, which 
has a continuous, rather than localized character. A better, although still schematic, 
picture is given in Fig. 2. As in the usual scattering processes reflection is due to the 
inhomogeneities of the potential, so here it is caused by the curvature. In fact, a region 
of spacetime containing a homogeneous gravitational field is flat, and it is known that 
waves are not diffused by a flat four dimensional background P, 0. 

As said above, this effect has been investigated for waves emitted by a compact 
source, by calculating explicitly the reflection and transmission coefficients, and has been 
found to be significant only in the induction zone and for long wavelengths, becoming 
negligible at large distances from the source and at late times, so that it is irrelevant 
for the problem of wave propagation in presence of astrophysical objects ||15|| . These 
general results can be qualitatively understood by remembering that scattering is non- 
negligible only in the regions where curvature is appreciable. Since the curvature of 
spacetimes associated to localized matter distributions drops off quickly as spatial infinity 
is approached, the region of scattering is localized near the source, as in Fig. 3. Moreover, 

1 We use units in which c = 1, and adopt the notations of Ref. M. The signature of the metric is +2. 



1 



since the "potential barrier" turns out not to be very high, most of the radiation escapes 
in the region where the curvature is negligible, and where it is no more backscattered. 
Reflection is appreciable only for long wavelengths. 

These considerations do not apply to spacetimes relevant as cosmological models, in 
which curvature is present over large scales and scattering could, in principle, take place 
everywhere (see Fig. 4). It is thus natural to ask whether in these cases the process can 
eventually "convey" a relevant part of radiation within the light cone; this constitutes 
the so-called cosmological tail problem [16]. Whereas qualitative criteria for the existence 



of tails have been already formulated and applied to special cases H-[[nj, no general 
quantitative characterization of the phenomenon, that could solve the problem above, 
seems to have been suggested up to now. A proposal of this kind, attempting to quantify 
the physical relevance of wave tails, will be described in a forthcoming article []17|] . In this 
paper we restrict ourselves to present (Sec. 2) a solution of Eq. ( |1 . 1|) in a form which can 
be useful for later references, and to point out (Sec. 3) the inapplicability of methods 
based on the calculation of reflection and transmission coefficients to investigate the 
occurrence of tails in a cosmological background. Sec. 4 contains some general remarks 
about the implications of this result. 

All the calculations will be performed in a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe, with metric?] 

ds 2 = a 2 (r]) [-drf + d X 2 + f(x) (d6 2 + sin 2 9 dip 2 )] , (1.2) 

where 

fsinhx x^(0,+oo) if K = -1 
X xe(0,+oo) if K = (1.3) 
sinx xe(0,7r) if K = +1 , 

and a(rj) is the scale factor. Although the general solution of Eq. ( |1 . 1| ) in a FLRW 
spacetime is well-known in an implicit form [18|, |2j, the explicit dependence on \ has 
been worked out, to the authors' knowledge, only in the case K — +1 []II| (the case 
K = is trivial). As far as the study of tails is concerned, however, this is the most 
intricate case, because the phenomenon, if present, is complicated by the possibility that 
radiation travels more than once through closed spatial sections. For this reason, and 
since it represents the only one not yet considered in the literature, we shall work out 
explicitly only the case K = — 1. 

2 rj is the so-called conformal time M. 
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2 General solution 

Equation Q) can be rewritten in the form 







-9 9 



ab 



8$ 



f.R$-m 2 $ = 0. 



w '—g dx a \ w a a <9a; 6 
In the FLRW background defined by Eqs. ( |1.2| ) and (|1.3|), Eq. ( |2.1| ) becomes 



(2.i; 



1 

P 



d 



dx \ dx 



+ 







1 <9 2 $ 



sin0 89 \ Sm9 89 ) + sin 2 9 d(p 2 



where now 



1 d (2^\ _ o 2^$_ a 2 m 2 $ = 

a 2 ar? V 8r\ J 



, a K 
# = 61 — + — 

a 6 a 1 



(2.2) 



(2.3) 



a dot denoting derivative with respect to rj. Separation of time and space variables, 



leads to 



and 



1 d ( 2 dT\ „ fc m 

^— a 2 — +£i?T + m 2 T+— T = 



9x \ dx) sm # \ ®0 j 



1 d 2 S 



+ kSf 2 = 



(2.4) 
(2.5) 
(2.6) 



sin 2 6> dtp 2 

where k is a separation constant. Equation (|2.5|) is in general difficult to solve, but it 
becomes very simple when considering a static universe |HJ [|, in which a =const. [] In 
this case the solutions are the usual exponentials, exp(±ia;ry), where 



k + a 2 m 2 + 6£K . 



Further separation of radial and angular coordinates in Eq. ( |2.6|) , 

S(x,0,p) = X(x)Y(9,p), 



(2.7) 



(2.8) 



The general solution of Eq. (|2.5|) is given in Ref. pfj also in the case a{rf) cx r] a (a / 1/2). 
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leads to the usual equation for spherical harmonics |21], |22 Y[ m (6, ip) and to 



d ( 2 <iX^ 



+ kf 2 X l = 1(1 + l)X t . 



(2.9) 



Setting ^i(x) = f(x)^i(x) an d using Eq. ( |1.3|) , Eq. ( p.9|) takes the form of a one- 
dimensional Schrodinger equation: 



dx 2 



+ 



2E 



1(1 + 1) 



where 



2E = k + K 



(2.10) 



(2-11) 



The explicit solution of Eq. ( |2.10| ) in the case K = +1 has been already given by 
Ford jl9| in terms of Gegenbauer functions. For K = 0, on the other hand, Eq. (|2.10|) 
becomes 



d 2 % 

dx 2 



+ 



2E- 



1(1 + 1) 

X 2 



#z = 



(2.12) 



which is familiar from the quantum mechanical description of a free particle in spherical 
coordinates (see, e.g., Ref. ||22|| ). Its general solution is a linear combination of 



^ 2 \X) = VX~H\ 



(1,2) , 
-l/2\ 



^Ex) 



(2.13) 



(1,2) 
1+1/2 



are the Hankel functions 21, 22 



where E > and H l 

In the case K = —1, Eq. Q2.10D allows the eigenvalue E to assume any positive value. 
It is convenient to define 



with p > 0; then 



dx 2 



+ 



2E = p 2 



2 W + l) 
P ~ 



sinh x 



% = 



Writing 

*Kx) = e**0,(x), 
Eq. ( |2.15|) transforms in the following equation for (frf 



d 2 <p, . dfa 1(1 + 1) 

+ 2 %p ; — tp. 



dx 2 



dx sinh 2 % 
4 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



It is now useful to define the new variable z e (—00, 0) as 

«=-(l-cothx) = r ^. 
In terms of z, Eq. (|2.17|) takes the form 

z(l-z)^ + (l-ip-2z)^ + l(l + l)<f> l = 
dz z dz 



(2.18) 



(2.19) 



which is immediately recognized as an hypergeometric equation |21| , P2j] with coefficients 



P = Z + l 
7 = 1 — ip . 

This equation admits two independent solutions [21, 2~2" ] , 

= F(a, (3, 7; z) = F{-1, 1 + 1,1 -ip; z) 

and 



(2.20) 



(2.21) 



t(2) 



r i-7(! _ z )i-<*-Pf(1 - p, 1 - a, 2 - 7; z) 



■1)V % F(-M + M + 



ip; z) 



(2.22) 



where F denotes the hypergeometric function. Correspondingly, Eq. ( |2.16| ) gives the 
two independent solutions of Eq. ( |2.15|) (an irrelevant factor (— l) ip has been dropped in 
getting Eq. from Eq. (ggg) ): 



*! 1) (x) = W(-MiM-ip(x)); 

*i (2) 00 = e-^FC-Z, Z + 1, 1 + ip; z{ X )) . 
The general solution of Eq. fl2.15| ) is thus 

*,(*) = 4 (p) *, (1) (x)+s,(p) ^ (2) (x) 



(2.23) 
(2.24) 

(2.25) 



with A[ and -B^ arbitrary complex functions of p. 

Notice that, since a is integer, the series defining F terminates, and therefore F 
reduces to a polynomial of degree I in z: 



F(-M + MTip;z) = £ 

n=0 



I + n \ (—l) n 



[1 =F ip)n 



(2.26) 
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where (A) n is defined for an arbitrary complex number A, by 

T(A + n) 



(A) n = A(A + l)-.-(A + n-l) 



r(A) 



(2.27) 



and (A) = 1. This result is useful if one wants to compare the exact solution ( |2.25| ) 
with the asymptotic one for \ —>■ 0. 



3 Asymptotic character of the solution 

In this section we argue that the concept of reflection and transmission coefficients is 
inapplicable to the study of tails in cosmological backgrounds. To this purpose, we 
shall restrict the discussion to the case K = —1 and a =const. For K = +1, due to 
the existence of closed spatial sections of spacetime, the transmitted radiation might 
be allowed to pass more than once through a given point of space, superposing to the 
fraction of radiation which is possibly reflected. For non-constant a(r)), the simple 
exponentials exp(±ipx) might not correspond to outgoing and ingoing waves. In both 
cases the treatment would suffer from unnecessary complications. 

The presence of tails is characterized by the fact that a pulse of radiation emitted at 
X = is partially backscattered. A stationary emission must therefore correspond to a 
solution of Eq. (|2.2|) which for x ~~ > contains incoming (reflected) as well as outgoing 
(emitted) radiation, whereas for x — * +o° only outgoing (transmitted) waves are present. 
Therefore, we must impose to the general solution ( |2.25| ) a boundary condition that 
corresponds to the absence of incoming waves at x ~~ * +oo. This is easily accomplished 
by noticing that for x ~~ *" +°° one has z — > and F — > 1. Hence, the asymptotic form 
of and #{ 2) is 

tfWft-^-KxOwe**, (3.1) 
M>z (2) (x — +oo) we"**, (3.2) 

which correspond, respectively, to outgoing and ingoing waves. The required boundary 
condition is therefore B\ = 0, leading to 

= Mp) e ip *F(-l, I + 1, 1 - ip; z( X )) (3.3) 

as the specific solution of the problem. 
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A calculation of the reflection and transmission coefficients along the same lines of 
that performed in Refs. |TT|-|[l^| would now require to expand *&i(x) as 



(3.4) 



where are independent solutions of Eq. ( |2.15|) corresponding to waves which are 
purely outgoing (+) and ingoing (— ) for x ~~ ¥ 0. Unfortunately, such solutions do 
not exist for / 7^ if we require purely ingoing/outgoing waves to have the form[] 
exp[— i{ujr] ± px)] as x 0- In fact, these simple exponentials are not solutions of 
the wave equation near the origin. This can be qualitatively understood by observing 
that in the neighbourhood of x — the centrifugal potential 1(1 + l)/x 2 varies extremely 
rapidly, preventing the modes exp(+ip X ) from occurring separately. As a consequence, 
one cannot select solutions of Eq. (|2.15|) which correspond to waves that near the origin 
have a purely ingoing or outgoing behaviour^]. 

Formally, this can be realized by characterizing these solutions as asymptotic eigen- 
functions of the operator P x = —id/dx, i.e., 



d^ 



(±) 



— 1 



+a + (higher powers of x) 



dx 



(3.5) 



with a > 0. Since for I 7^ 



H h P x =O( X - s )^0 



(3.6) 



where 



1 d 2 



i(i + i) 

2 dx 2 ^ 2 sinh 2 x 



it follows that there are no common eigenstates of the operators Hi and P x , i.e., 
are no solutions of Eq. ( j2.15 ) which satisfy the asymptotic condition ( |3.5| ). 



(3.7) 



there 



4 The following results are conditioned by the use of this characterization, which appears rather 
restrictive. However, other definitions do not seem sharp enough in selecting waves which can be 
considered as purely ingoing or outgoing. The formulation of a more general characterization would 
nevertheless constitute an important progress as far as the calculation of reflection and transmission 
coefficients is concerned. 

5 It is easy to check that the same conclusion holds for a generic solution of the wave equation that 
contains components with / =/= 0. An analogy with the behaviour of particles can be drawn starting from 
the similarity of Eq. (2.15) with a Schrodinger equation in a central potential. The intuitive explanation 
is that particles with nonzero angular momentum cannot "enter in" or "exit from" the centre \ = 0. 
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This conclusion can be also checked explicitly as follows. Looking for solutions of 
Eq. ( p. 15 ) of the form 

Vi(x)=X m + 0( X m+1 ), (3.8) 

one finds immediately that either m = — I or m = l+l. The general solution of Eq. ( |2.15| ) 
can therefore be written as 

Mx) = C l f l (x)+D l g l ( X ), (3.9) 

with Ci and Di arbitrary complex numbers and fi, gi two particular solutions which 
behave, for X 0, as 

f l ( x ) = l + 0(x 1 - 1 ), (3.10) 
X 

and 

9i(x) = X W + 0( X l+2 ) • (3.11) 
The coefficients Ci and Di for any particular solution can be obtained by comparing 
Eqs. (|3.9|) - (|3.11|) with the expansion of ^i{x) f° r X ~~ * 0. In particular, for one can 



write 

M ± \x) = Ci ±) f l (x) + Dl ±) g l (x), (3.12) 
where the coefficients and Di should be chosen in such a way as to guarantee 
the prescribed character (purely outgoing or ingoing) of the solutions for x 0. 
Substituting Eq. ( E>.12| ) into ( |3.5| ) and requiring the coefficients of the leading powers of 
X to vanish, one obtains the trivial result C\ = d\ = 0. 

The non-existence of solutions of Eq. (|2.15| ) corresponding to waves which are purely 
ingoing or outgoing in the region x — > 0, prohibits one to define reflection and trans- 
mission coefficients as in the usual treatments of scattering problems. It seems that no 
method has been developed in the literature to deal with similar situations. The case 
considered here is quite different from the simpler ones arising when studying diffusion 
of waves by Schwarzschild black holes, where a coordinate transformation can be found 
for which the corresponding one dimensional Schrodinger problem involves a finite, lo- 
calized potential barrier, and nontrivial asymptotic eigenfunctions of momentum can be 
found PH-P . 



4 General remarks 

The most appropriate treatment of the tail problem would be the explicit determination 
of the reflection coefficient describing backscattering by the cosmological curvature, and 
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characterizing quantitatively the fraction of radiation which does not propagate along 
the light cone. The impossibility of carrying on this approach leads one to look for 
alternative ways of studying the phenomenon. A straightforward idea could be to pursue 
the formal analogy between the tail problem and quantum scattering, by observing that 
Eq (|2.15|) can be thought of as the radial part of a stationary Schrodinger equation 



V (X) = — — ~ 72 ■ ( 4 - 2 ) 



describing a particle with Hamiltonian operator 

h = Ip 2 + 1 -v(x)1 2 , (4.1) 

where p and L are, respectively, the linear and angular momentum operators in a ficti- 
tious three-dimensional Euclidean space in which x plays the role of an ordinary radial 
coordinate, *&i(x)/x is the radial part of the l-th component of the Schrodinger wave 
function, and 

1 1 
sinh 2 x X 2 

One possibility would then be to compute the cross section for this quantum mechanical 
system. Although such a calculation could be carried on, it cannot however be regarded 
as a satisfactory solution to our specific physical problem, which concerns waves emitted 
from x — rather than incoming from infinity. A satisfactory study of the subject seems 
thus to require a radically different approach. In Ref . |T7| , a technique will be presented 
which characterizes tails in terms of the ratio between their energy content and the total 
energy of the field. 

The impossibility of defining reflection and transmission coefficients is not the only 
reason (although obviously a compelling one) to look for a different characterization of 
tails. Even if such coefficients could be defined, in fact, they would not account in a 
reliable way for the fraction of radiation propagating inside the light cone. This can 
be realized by noticing that a crucial part of the information about the character of 
wave propagation is exclusively contained in Eq. (|2.5| ), whereas the calculation of the 
reflection and transmission coefficients involves only Eq. fl2.9| ). For example, as far as 
Eq. (|2.9|) (and, hence, the reflection and transmission coefficients) is concerned, all the 
values of £ and m are equivalent, which is certainly not the case since, e.g., the choice 
m = 0, £=1/6 leads to no tails in every conformally flat spacetime ||. A way to take 
into account the time dependence as well is to perform the separation 

^x,e,<p) = ^r^-y(e,v>) (4.3) 
fix) 
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in Eq. (|2^), getting 



d 2 ipi d 2 i>i 



1 



dx 2 dr] 2 
Under the coordinate transformations 



6£ + a m H * — 

sinh x 



(4.4) 



w = | fa - x) 

\{v + x) 

u = tanh u 
v = tanh v , 



r/ 

X z 



- v + u 
v — u , 



Eq. (|4.4D becomes 



c?x 2 dfj 



1 1 2 



l-6^ + a 2 m 2 Z(Z + 1 



M 2 )(l 



+ 



X 2 



(4.5) 



(4.6) 



(4.7) 



(4- 



If 



1 - 6 f + a 2 m 2 



, (4.9) 

Eq. (|4.8|) , and consequently Eq. ( p.l|) , is clearly tail-free, since it corresponds formally to 
wave propagation in flat spacetime[]. This agrees with previous literature || ||, and can 
be related to the known fact that sinh -2 x is a so-called reflectionless potential |23|, 24 . 
However, both this argument and previous works are not relevant when Eq. (|4.9|) does 
not hold - for example if m = £ = 0. 

In spite of the deficiencies of the treatment, one feature of the phenomenon emerges 
nevertheless clearly: The tail problem for radiation in FLRW spacetimes regards in 
general only waves whose wavelengths are at least of order a(i]). This can be seen 
explicitly for the case of a static universe with K — —1, in which Eqs. (|2.7|) , ( 2.11 ) and 
( |2.14j ) give, for the case m 



0. 
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p 2 + l 



6£ 



(4.10) 



that can be regarded as the dispersion relation for the waves. From Eq. Q4.10] ) it follows 
that, except in the very special case £ = 1/6, only components with frequency and wave 



3 We thank an anonimous referee for pointing out this consideration. 
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number[] smaller than (or comparable to) a" 1 are subject to diffusion. Since tails would 
then regard only a very small, extreme part of the spectrum of the radiation content of 
the universe, and since the variability of these waves on scales of order a is too slow to 
be detected, it appears that the effect is hardly observable. However, no conclusion can 
be drawn without a careful investigation; in fact radiation of very long wavelengths can 
possibly lead to physically relevant effects (see, e.g., Ref. |25|). 

Scalar fields satisfying a wave equation are considered in the inflationary models of 



the early universe (see, e.g., |26| and references therein). The tail problem for the inflaton 
regards length scales relevant for cosmology, so it could possibly have some importance 
in problems connected to the physics of the early universe. This is however beyond the 
scope of the present work. 
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